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by 

K.  W.  H.  Stevens 

Cruft  Laboratory,  Harvard  University 
Cambridge,  Massachusetts 

Abstract 

Previous  attempts  to  explain  the  anisotropy  of  the  quadrupole 
resonance  in  solid  iodine  are  critically  reviewed  and  an  attempt  is  made 
to  set  up  a valence  bond  description  of  the  iodine  molecule  which  con- 
tains the  possibility  that  additional  weak  bonds  are  formed  in  crystalline 
iodine.  This  is  primarily  accomplished  by  introducing  d-hybridization . 
The  observed  bond  angles  are  then  used  to  estimate  some  of  the  param- 
etersinthe  orbits.  Use  is  made  of  knowledge  about  nuclear  quadrupole 
interaction  in  solid  in  the  iodine  atom  to  determine  further  pa- 

rameters, and  in  this  way  a picture  of  the  iodine  molecule  and  its  be- 
havior in  a crysteil  is  obtained.  An  interesting  possibility  which  then 
arises  is  that  the  lobes  of  the  wave -function  of  different  iodine  molecules 
do  not  overlap  very  strongly  because  a maximum  of  one  molecule  seems 
to  be  directed  towards  a minimum  of  another.  This  does  not  correspond 
to  the  usual  ideas  about  bonds,  but  as  the  theory  passes  smoothly  from 
the  notion  of  maximum  overlap  to  this  new  situation,  no  discussion  is 
given  about  how  it  should  be  described.  It  is  also  seen  that  the  observa- 
tion of  the  quadrupole  resonance  in  single  crystals  of  iodine  (if  this  be- 
comes possible)  would  provide  very  useful  information  and  checks  on  the 
theory. 
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The  observation  of  pure  nuclear  quadrupole  resonances  in  solid  iodine 
has  been  reported  by  a number  of  authors  [1],  who  find  that  there  are  con- 
siderable departures  from  axial  symmetry  in  the  charge  distributions  near 
the  iodine  nuclei.  At  first  sight  this  is  rather  surprising,  for  iodine  forms 
a molecular  crystal  of  molecules,  and  ^ne  would  expect  that  the  charge 
distributions  would  show  cyliiidrical  symmetry  about  the  I-I  bond,  as  is 
certainly  true  for  a free  molecule.  An  explanation  of  the  origin  of  this 
asymmetry  has  been  put  forward  by  Townes  and  Dailey  [2],  who  point  out 
that  the  I-I  bond  length  in  the  solid  is  greater  than  in  the  free  molecule,  a 
fact  which  can  be  understood  if  additional  weak  bonds  are  formedinthe  crystal, 
thus  weakening  the  main  bond.  A study  of  the  crystal  structure,  see  Figs,  1 
and  2,  shows  that  the  iodine  molecviles  lie  in  planes  which  are  fairly  well 
spaced,  and  that  in  the  planes  there  are  two  directions  for  the  I-I  bond,  the 
different  sorts  of  molecules  lying  in  chains.  In  Fig.  2 we  show  an  enlarged 
picture  of  an  iodine  molecule  with  the  positions  of  a pair  of  nearest  neighbor- 
ing iodines,  and  if,  following  Townes  and  Dailey,  we  assume  that  the  additional 
bonds  are  formed  in  these  directions,  it  seems  fairly  clear  that  the  charge 
distribution  near  an  iodine  nucleus  will  not  have  an  axis  of  symmetry.  The 
rest  of  the  treatment  of  this  problem  by  Townes  and  Dailey  is  more  difficult 
to  understand.  They  assume  that  the  charge  distribution  around  each  of  the 
three  bonds  made  by  any  iodine  has  cylindrical  symmetry  about  the  bond,  so 
that  only  one  parameter  is  needed  to  describe  the  quadrupole  interaction  with 
the  charge  in  each  bond  (if  z is  a bond  direction  d^dx^  = 3^/8y^=  8^/dz^), 

and  then,  because  the  two  weak  bonds  are  of  equal  length,  they  weight  them 

* On  leave  from  Nottingham  University,  England 
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each  with  a total  charge  a and  the  main  bond  with  l-2a.  In  this  way  the  total 
quadrupole  interaction  of  an  iodine  nucleus  is  described  by  two  parameters, 
a and  8^V/8z^;  a is  found  from  the  anisotropy  and  8^V8z^  from  the  magnitude 
of  the  interaction.  The  final  step  is  thea  to  show  that  8^/8z^  found  in  this 
way  agrees  with  the  field  gradient  expected  for  a free  iodina  molecule.  Our 
main  criticism  of  this  procedure  is  that  there  is  no  reason  to  suppose  that  be- 
cause the  two  additional  bonds  are  equal  in  length  they  must  have  equalweight, 
nor  is  there  any  a priori  reason  to  assume  that,  apart  from  the  different  charges, 
8^/8z^  is  the  same  for  the  weak  and  strong  bonds.  This  latter  assumption  be- 
comes plausible  if  one  assumes  that  the  bonds  are  made  of  hybridized  s and  p 
orbitals,  but  there  is  then  the  difficulty  that  three  bonds  in  a plane  cannot  be 
made  from  s and  p alone.  The  first  criticism  can  be  illustrated  very  simply 
(see  Fig.  3)  where  the  dotted  lines  now  represent  the  charge  densities  in  the 
additional  bonds  in  a schematic  way.  The  bonds  will  clearly  have  equal  lengths 
but  they  do  not  set  up  equal  field  gradients  at  a given  iodine:  we  believe  that 
this  sort  of  behavior  is  more  likely  to  occur  than  equal  gradients. 

With  these  difficulties  about  the  validity  of  Townes  and  Dailey's  model 
we  decided  to  try  to  set  up  a more  specific  model,  based  on  the  Paulingideas 
of  directed  valence.  What  is  needed  is  a reasonable  description  of  a free 
iodine  molecule,  then  a study  of  how  such  molecules  might  come  together  to 
form  a crystal  and  finally  a calculation  of  the  nuclear  quadrupole  interaction. 
This  ambitious  programme  proves  to  be  too  difficult  to  carry  through  com- 
pletely, but  quite  a lot  of  progress  can  be  made  and  a number  of  useful  results 
obtained. 

Since  this  work  was  completed  another  theory  about  the  origin  of  the 
anisotropy  has  been  presented  by  Robinson,  Dehmelt  and  Gordy  [3]  . This 
differs  significantly  from  the  previous  theory  and  from  the  one  we  propose 
to  give.  While  we  do  not  know  which  theory  is  nearest  the  true  solution,  we 
feel  that  ours  has  at  least  a more  fundamental  approach  than  the  other  two 
and  that  something  of  this  sort  is  required,  rather  than  further  semi- 
empirical  speculations . 
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FIG.  I THE  ARRANGEMENT  OF  IODINE  MOLECULES  IN  A PLANE 
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FIG.  3 AN  EXAMPLE  WHICH  SHOWS  SCHEMATICALLY  HOW 
BOND  DISTANCES  CAN  BE  EQUAL  BUT  GIVE  RISE 
TO  DIFFERENT  CHARGE  DISTRIBUTION  AT  A GIVEN 
NUCLEUS 
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The  Free  Molecule 


Each  iodine  atom  has  five  5p  electrons  outside  closed  shells  and  we 
might  expect  that  in  the  molecule  four  of  these  pair  off  in  p^  orbitals  and 
the  remaining  one  goes  into  a orbit  and  pairs  up  with  a similar  electron 
from  the  other  iodine.  This  is  too  simple  because  hybridization  is  expected 
and  to  allow  for  this  we  shall  generalize  the  above  orbits  by  expressing  them 
in  a series  of  spherical  harmonics.  Thus,  if  the  positive  direction  of  z is 
taken  towards  the  other  iodine,  we  have: 


2 electrons  in 


2 electrons  in 


1 electron  in 


where  allowance  is  made  for  s,  p,  d hybridization,  [yz^etc.,  denotes  the 
angular  variation  of  a normalized  wave-function,  and  N2  are  normalizing 
factors  and  a,  a and  y describe  the  relative  admixtures  of  orbitals  of  dif- 
ferent symmetries.  The  normalization  is  such  that  N2  will  probably  be  near 
unity  (not  exactly  because  this  orbit  overlaps  strongly  with  the  similar  one  on 
the  other  iodine),  while  N ^ will  depart  from  unity  according  ; s a is  large  or 
small  > 

The  bonding  orbit  (i.e.,  the  orbit  with  only  one  electron)  is  expected 
to  have  a strong  maximum  in  the  positive  direction  of  z,  but  this  is  not  the 
only  direction  of  maximum  charge  concentration,  particularly  with  d-hybridiza- 
tion.  In  fact,  if  6 is  measured  from  the  positive  direction  of  z the  angular 
variation  of  the  amplitude  of  the  wave-function  at  some  radius  r is  proportioned, 
to: 

a/s’  2 

u cos  aV3  cos  0 + v sin  a ( 3 cos  6-1) 

where  u and  v denote  the  values  of  the  radial  parts  of  the  p and  d orbitals 
respectively  (u  and  v are  usually  taken  to  be  equal).  ■ This  function  hasturn- 
ing  values  at  9 = 0,  180°  and  6 - cos  < — cot  a | . To  make  6 = 0 . a 

0,1  J 
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strong  maximxim  cot  a must  be  positive,  so  that  0^  lies  between  90°  and  180°. 
Thus  as  a result  of  using  d-hybridization,  the  bonding  orbit  has  subsidiary 
turning  values  on  the  cone  of  directions  corresponding  to  0^  as  well  as  at 
0 = 180°.  It  is  convenient  to  have  a notation  for  the  directions  of  these  turn- 
ing values,  so  we  shall  refer  to  them  by  the  angles  they  make,  i.e.,  0°,  180°, 

and  6 directions, 
o 

Crystalline  Iodine 


Suppose  two  iodine  molecules  are  coming  together;  then  we  may  ask 
how  they  will  arrange  themselves.  If  we  neglect  for  the  moment  the  disturb- 
ances of  their  charge  clouds  by  their  mutxial  interactions,  it  seems  not  un- 
likely that  they  either  have  their  0^  directions  pointing  towards  one  another, 
or  that  the  0^  of  one  will  point  towards  a 180°  direction  of  the  other.  Referring 
to  the  crystal  structure  (Figs.  1 and  2)  it  is  seen  that  one  of  the  subsidiary 
bonds  is  only  10°  away  from  the  180°  direction  and  that  the  other  lies  between 
90°  and  180°.  Thus  a s emi -empirical  understanding  of  why  iodine  has  the 
structure  it  does  is  obtained  if  we  assume  that  the  molecules  arrange  them- 
selves so  that  the  180°  direction  of  one  is  directed  towards  the  0 direction 

o 

of  another.  Thishypothesis  will  be  further  strengthened  if  we  can  show  that 
the  mutual  interactions  of  the  molecules  reduce  the  strength  of  the  main 
I-I  bond  and  move  the  180°  direction  in  the  right  direction.  To  consider 
this  we  suppose  that  the  interactions  cause  mixing  of  our  original  three  orbits, 
so  that  in  the  crystal  the  electrons  are  distributed  in  orbits  as  follows: 


2 in 


2 in 


1 in 


cos 


Nl 

cosg 


|y>  + alyz>j 

|x>  + a|  xz>]  +-^ 


cos  y ( cos  o I + sin  a | 3z 
+ sin  V ! 


N. 


3 Z 1 

cos  "Y  (cos  a 1 ^ + sin  o j 3z^  - r + sin  y | 
-T^[|x>+a|xz>j 


where  the  plane  of  the  iodine  molecules  is  taken  as  xz,  so  the  orbit  which  is 
primarily  perpendicular  to  this  is  unaltered.  (We  are,  for  the  present. 
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neglecting  the  3-dimensjQnal  nature  of  the  problem  and  considering  only  a 
plane  of  iodine  molecules.).  If  we  assume  that  p is  small,  the  turning  points 
of  the  modified  bonding  orbital,  in  the  xz  plane,  are  along  9 = (|4,  tt  + 
i 9^  + X where 


-tang  ^2  (u  + a ) 

cos  y (ucos  a + V 15  v sin  a) 

- tan  B ^2  (u  - a vV?) 

cosy  NT  , *rTc"  \ 

* 1 (u  cos  a -V15  V sina) 


andX 


* N-  AV3ucos9  +av^Jr5cos29  I 

- ~l3.n  p 2 L ’ o ^ ^ 


cos  Y N 


1 [VJu  cos  Q cos  0^  + sin  a cos  29^j 


With  some  rearrangement,  and  using  the  definition  of  9^,  these  become: 

. ^ CO8  0 (1+A) 

j _ - tan  6 o'  ' 

X (cose  - 1) 

' o 


‘I' 


. - cos  0 (1-A) 

-tanp o' 

X (cos  6 + 1) 

' o 


X 


- tanp 


•ol' 


cos  0 (cos  9 + A cos  20 


1 - cos  0 


[cos  9 + A cos  29  { 
I o o3 

(Y-co.-e-yn'-’Ar 


where 


X = cos  a cos  y 


N. 


and 


A = a1[5  X 

' u 


P is  positive  if  the  bond  along  + 9^  is  strengthened,  but  whether  the  180° 
bond  is  bent  in  the  required  direction  depends  on  whether  A is  less  than, 
or. greater  than,  unity.  Alternatively  v.'e  can  identify  41  with  its  experi- 
mental Vcilue,  10°53’  and  0^  + X similarly  with  105°25  * and  use  the  above 
equations  to  determine  0^  and  tanp/Xas  functions  of  A.  The  results  of  this 
are  given  in  Table  1. 
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Table  1 


Values  of  parameters  needed  to  fit  the  bond  angles 


A 

0 

.2 

.4 

.6 

9 

o 

108° 

110° 

113.7° 

120° 

tan  3 

.425 

.457 

.471 

.475 

The  fact  that  is  different  from  zero  means  that  the  main  iodine  bond  is  not 
directed  straight  towards  the  other  iodine,  which  suggests  that  the  main  I-I 
bond  is  weakened  in  the  crystal,  as  required.  Further,  to  have  a stable 
arrangement,  the  main  bonds  from  the  two  iodines  in  a molecule  must  be 
deflected  on  opposite  sides  of  the  I-I  axis,  which  means  that  the  midpoint 
of  a molecule  is  a center  of  symmetry. 

One  point  which  has  been  glossed  over  and  which  must  now  be  considered 
is  the  use  of  the  turning  values  of  the  angular  part  of  the  wave -function  to 
determine  bond  directions.  In  so  far  as  the  Pa;iling-directed-bond  method  is 
valid,  it  would  seem  reasonable  to  associate  bond  directions  with  maxima, 
but  in  the  above  reasoning  we  have  not  shown  that  our  various  directions, 

0,  0^  and  180°  are  really  directions  in  which  the  electron  distribution  is 
large.  In  particular  it  seems  not  unlikely  that  the  actual  orbits  will  be  such 
that  the  charge  in  the  180  direction  is  very  small,  for  presumably  the  bulk 
of  the  charge  is  located  in  0 = 0.  If,  then,  we  assume  that  the  charge  concen- 
tration is  a minimum  in  0 = 180  we  may  picture  a free  iodine  molecule  as 
having  a lot  of  charge  along  the  main  axis  and  with  subsidiary  cones  of  charge 
at  about  105°  with  this  axis.  The  crystal  is  then  built  up  by  having  a point  on 
the  cone  of  one  molecule  on  the  axis  of  the  cone  of  another.  In  this  way  the 
molecxiles  pack  fairly  closely  without  their  charges  overlapping  too  strongly, 

1. e.,  without  the  formation  of  extra  bonds.  It  also  seems  not  unreasonable 
to  suppose  that  this  would  lead  to  planes  of  fairly  closely  interlocked  iodine 
molecules,  with  weaker  forces  between  the  planes. 

The  above  ideas  are  based  on  physical  arguments  about  charge  densi- 
ties, and  it  is  seen  that  the  introduction  of  d -hybridization  is  mainly  responsible 
for  the  results  obtained.  It  is  clearly  going  to  be  a formidable  task  to  try  and 
substantiate  these  ideas  by  actual  calculation  of  orbits,  so  we  propose  instead  to 
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investigate  the  results  of  the  nuclear  quadrupole  resonance  to  see  what  in- 
fornnation  this  can  give. 

The  Nuclear  Quadrupole  Interaction 

If  (x,y,z)  is  the  position  of  an  electron  with  respect  to  a nucleus  at  the 
origin,  the  contribution  which  it  ma.kes  to  the  electric  quadrupole  interaction 
is: 


e^Q 


^(I  1 +I  1 )+^^(l  1+11) 
5 yz  zy  5 zx  xz 
r ^ ^ r 


To  obtain  the  full  interaction  this  must  be  summed  over  all  electrons  and 
averaged  over  their  orbits.  Using 


r \r  / d 


✓ ^ zx  2 2 x_  1 / 1 \ 

l3z  -r  >—H) 

T 7yT\r/d 


and  other  quantities  derived  from  them,  the  interaction  in  the  present  case 
become  s : 


e^Q 


21(21-1) 


a,, I ^ + a,,l  ^ + a-_l  ^ + a,,(I  1 +1  I ) 
11  x 22  y 33  z 13  z x x z 


where 


■11 


4 1 1,2,  2 cos^Y  /I  2 ^1,2, 

z(+^-=\a  ) ^(_cos  a+-\sin  a) 

Nj2  57  N2^  5 7 

+ sin^  — ^ (^  + ^X.a^)  + ^ — 3£.(-^cos^a  -^X,sin^  a) 


t 
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r 


4,1  1,2.  2co«  V.  1 2 . 1 , . 2 


22 


-(  1 \ 


N, 


N 


( ^ cos  a +=-  Xsin  a) 


2 

+ sin  p 


1 , 2 


IN 


2 , 1 2,2..  2 cos  Y , 1 2 i , . i,  .1 

( - ^ - y Xa  ) + .2  ' (-^cos  o-yXsin  o)| 


N. 


‘3,  A) 


N 


4 , 2.2.2.  2 cos  Y/  2 2 2 . • 2 . 

( - ^ + y Xa  J 1”"*"'  " s’  “ "T  a ^ 


N. 


^ . 2-  1 -2  , 11-2.^  2cos  Y,  2 2 ^ 2 . . 2 . 

+ sin  6 < «(  - ^+  — Xa  } + ■ ^ cos  o + Xsin  of 

'^257  25  7 

^1  ^2 


13' 


- f _L_)  2 cos  p sinp  cos  Y 


Ni«2 


cos  a 


sin  a 


where 


These  formulae  are  too  complicated  to  be  handled  easily,  so  as  ( — 5) 

1-3  ^ 

probably  varies  as  t.  for  electrons  with  the  same  principal  quantum 

number,  it  seems  reasonable  to  approximate  by  taking  \ as  zero.  Con 

siderable  simplification  follows  if  ratios  are  then  used,  giving: 


ail  1 - 7^  - sin^p  (2+X^f 
^33  -2+  X^+  j sin^p  (1  + 2X^f 

^22  1 - i sin^  (1-X^f 

^33  -2+  X^+  isin^p  (l  + 2X^f 

111  = 1 cosB  sinp  X 

^33  ^ -2+  X^+  j sin^p  (l  + 2X'^f 


L 
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These  relations  are  not  independent  since  ^ ZZ  ^ ^33  ~ ^ 

investigation  of  the  quadrupole  interaction  would  determine  p and  X,  and 
then  (see  Table  1)  A and  0^  which  are  properties  of  the  free  molecule. 

There  are,  however,  a number  of  experimental  difficulties  in  this,  so  that 
present  information  extends  only  to  the  magnitude  of  and  the  anisotropy 
I Ajj^  - capital  A's  mean  that  the  quadrupole  interaction 

is  referred  to  its  principal  axes  (i  . e . , A ^ ^ ~ assume  that  a^^is 

small  so  that  the  original  axes  are  close  to  the  principal  axes,  the  anisotropy 
is  approximately  | 3-^]^  " ^2zl/^33 

SAn  p 

I -2  + j sin^p  ( 1 + 2X^)| 

To  get  the  bond  angles  right,  we  see  (Table  1)  that  tan  p /x  is  probably  of 
order  .45,  so  taking  tan^p  = .18X^  the  observed  value  of  .15  for  the 
anisotropy  gives  X^^x/  .75  and  tanp  r\j  .37.  Thus,  although  the  theory  con- 
tains a good  many  parameters,  some  numerical  values  can  be  obtained 
with  some  precision.  Continuing,  with  X^  = .75,  N cos  a cos  y = . 866, 
and  using  cos  0^  = -l/^fTS  u/v  cot  a and  setting  u = v,  0^  = 110°,  we  obtain 
a /X/  37  andNj/N^cos  y/xyl.l.  If  the  nonbonding  orbits  contain  no  d-wave 
functions,  is  probably  less  than  and  this  last  relation  will  not  be 
satisfied.  It  therefore  seems  that,  with  our  model,  it  is  necessary  to  have 
enough  d-character  in  the  nonbonding  orbits  to  make  and  also  we 

may  expect  y-^O.  With  a/xy37°,  yrvyO  the  bonding  orbit  contains  very 
little  s-character  and  approximately  equal  amounts  of  p and  d,  which  means 
that  0 = 180°  is  a direction  in  which  there  is  very  little  charge.  In  other 

words,  the  p and  d admixtures  are  such  that  the  p approximately  cancels 

o 

the  d in  the  0 = 1 SO  direction.  This  situation  was  anticipated  in  the  dis- 
cussion earlier  in  the  paper.  A more  surprising  result  is  that  the  tr  -orbitals, 
i.e.,  those  which  contain  pairs  of  electrons  moving  perpendicular  to  the 
axis  , have  /x;l0%  d-character  (amplitude  of  admixture, from  Nj),  but  on 
examination  this  is  seen  tc  make  the  electrons  concentrate  towards  the 
middle  of  the  molecule,  which  does  not  seem  unreasonable. 

The  above  numerical  estimates  were  obtained  by  taking  X.  = 0 and 
neglecting  If  we  new  use  the  above  numerical  values  and  do  not  make 
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these  assumptions  it  is  found  that  the  theoretical  anisotropy  might  be  re- 
duced to  . 14,  depending  on  the  value  of  The  observed  value  can  be 
restored  by  increasing  X slightly,  a change  which  will  not  affect  the  main 
conclusions . 

Finally  , it  is  of  interest  to  consider  the  actual  magnitude  of  the  quadru- 
pole  interaction.  Our  equations  give: 


^^33^  crystal 
^ ^ 3'3^free  mo'^ 


YN-/ 


;2  + X^  + j sin^P  (1+2X^) 
— — 2 


01^1 


r\J  . 87  6 


The  measured  value  of  the  quadrupole  coupling  in  the  solid  is  2140  Mc/s 
so  this  ratio  predicts  that  the  free  iodine  value  should  be  2440  Mc/s  . 
The  other  thing  which  can  be  done  is  to  compare  (^33)j.yygtai  with  the 
corresponding  value  for  a free  iodine  atom,  since  a free  atom  will  oe 
described  by  letting  a = P = Y = a = 0 and  Nj  = N2  = 1. 

2 12  2 

(a--)  -2  + X'^+~sin'^P(l+2X'^)  , 

33  crystal  2 ' _1.Z3 

"I  ■ '■  ""T‘  " * ’ I -I  I.  y II  I I II  — S 

'^33'free  atom  -2N,  N- 


1 


1 


(a,_),  can  be  obtained  from  the  analysis  of  the  hvperfine 

33  free  atom  2 1 23  2292 

structure  of  free  iodine  by  Jaccarino  [4],  whence  = 1. 32. 

The  departure  of  from  unity  gives  a measure  of  the  amount  of  d- 

character  in  the  nonbonding  orbit,  and  this  agrees  well  with  our  earlier 

estimate  of  this. 


Conclusion 

Starting  with  elementary  ideas  about  the  nature  of  chemical  bonds,  it 
has  proved  possible  to  set  up  a plausible  model  of  a free  iodine  molecule 
which  enables  one  to  understand  how  iodine  molecules  may  come  together 
to  form  a solid  having  the  observed  structure.  The  theory  is  then  tested 
by  examining  the  predictions  which  it  makes  about  the  nuclear  quadrupole 
interaction.  A number  of  approximations  are  made  in  this,  but  it  is  found 
that  the  anisotropy  and  magnitude  of  the  quadrupole  interaction  can  be 
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correlated  in  a satisfactory  way  with  the  crystal  structure  and  wave 
functions,  thus  helping  to  confirm  the  initial  assumptions.  The  final 
picture  of  the  origin  of  the  anisotropic  quadrupole  interactions  then 
differs  significantly  from  the  models  of  Townes  and  Dailey  and  of 
Robinson,  Dehmelt  and  Gordy. 
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